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Using the spin-wave approximation elementary excitations of a semi-infinite two- 
dimensional -S = 5 Heisenberg antiferromagnet are considered. The spectrum consists 
of bulk modes - standing spin waves and a quasi-one-dimensional mode of boundary 
spin waves. These latter excitations eject bulk modes from two boundary rows of sites, 
thereby dividing the antiferromagnet into two regions with diff'erent dominant excita- 
tions. As a result absolute values of nearest-neighbor spin correlations on the edge exceed 
the bulk value. 
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1. Introduction 

It is well known that at certain conditions defects of crystal structure can generate 
states localiz ed ill the defect region, while the magnitude of bulk states is suppressed 
in this region.SEI This leads to the situation in which the defect neighborhood and 
the rest of the crystal constitute two systems with different dominant excitations. 
A crystal surface can be also considered as a defect.^ As applied to the surface, 
the mentioned situation leads to the appearance of a near-boundary region whose 
properties differ from the bulk properties. 

The influence of boundaries on the spectrum a nd ob servables of the quantum 
Heisenberg antiferromagnet has been studied in 

twcPSH (2D) and threeP (3D) di- 
mensions. In particular it was shown that absolute values of the nearest-neighbor 
spin correlations near the boundary exceed the bulk value. This result was in- 
terpreted as a manifestation of increased valence-bond-solid correlations near the 
edge.EEl 

In this Letter, we propose another interpretation of the increased spin corre- 
lations near the edge. We relate their appearance to the peculiar spectrum of the 
semi- infinite antiferromagnet. The spectrum involves bulk modes - standing spin 
waves and a quasi-one-dimensional mode of boundary spin waves. These latter ex- 
citations are observed in the two boundary rows of sites, and they eject the bulk 
excitations from this region. Thus the antiferromagnet appears to be divided into 
two regions with different dominant excitations. As known, nearest-neighbor spin 
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correlations of the one-dimensional (ID) antiferromagnet are larger than in the 
2D case. As a consequence the quasi-lD near-boundary mode prod uces l arger spin 
correlations than the 2D bulk modes, which explains the observecPEEl increased 
near-boundary correlations. Similar interpretation can be applied to the 3D case.^ 
As will be seen below, the description of perturbations introduced by the edge into 
the magnon sp ectru m is in many respects similar to the problem of a local de- 
fect in a 

crystal.lIEl Thus seemingly different imperfections of crystal structure are 
described in the framework of the same approach. 

2. Model and its Elementary Excitations 

We suppose that the edge is located along one of the crystallographic axes and 
choose coordinates so that the antiferromagnet is situated in the half-space Ix > 
and described by the Hamiltonian 

where sites of a 2D square lattice are labeled by the two coordinates Ix and ly , the 
lattice spacing is set as the unit of length, and 8;^^;^ is the spin-^ operator. 

For the temperature T = the antiferromagnet has the long-range order and 
its elementary excitations can be described in the spin-wave approximation, 

5f = e'™Q-6^), S^^P,H, + P^bl (2) 

where 11 — (tt, tt), P^ = i (l ± e*'^'), and 1 = {Ix, ly). The spin-wave operators bi 
and &f satisfy the Boson commutation relations. Substituting Eq. ([2|) into Hamil- 
tonian ([1]) , dropping terms containing more than two spin- wave operators and con- 
stant terms, we obtain 



H^JJ2T. 2 - i'^'^") ^lyi.^i^ + cos(fcy) {bkyi.b-ky,i. + bli/^ky,i. 

+ {bk^lj-ky,l^ + l + &I^,^6Lfc^,;^ + i) 



(3) 



Here the translational invariance of the Hamiltonian along the Y axis was taken 
into account and the Fourier transformation bk^i^ = N^^/"^ e^'^^^^bi^i^ was used 
with N the number of sites in the Y direction and ky the ID wave vector. 

To investigate the spectrum of elementary excitations we introduce the two- 
component operator 

and define the matrix retarded Green's function 

D{kytlxl'x) = -ie{t) ( [Pk^i^ (t), Bl,^] ) . (4) 




I 



I 
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In Eq. (gl), Bkyi^it) = e'-^*i3fe„/,e~*^* with the Hamiltonian determined by Eq. ^ 
and the angular brackets denote the statistical averaging. 

To calculate Green's function ([4]) we use the equation of motion, 



i-D{kytlJ'J = 3{t)6i^i,j3 + J 



where 



2(1- -^Si^o ) fa + cos(fcy)fi 



b{kytlX) 



J. 

Ti 

2 ^ 



b{kyt, I, + 1, Q + b{kyt, I, - 1, Q 



(5) 



ri = 



1 

-1 



T3 = 



1 

-1 



The solution of Eq. ^ can be written in the form 



D{kytlj'J = b^'^KkytlX) - - I dt'b'^''\ky,t-t',l,0)b{kyt'OQ, (6) 



where Z)*^"^ {kytlxl'x) is Green's function corresponding to Hamiltonian ^ with- 
out the term proportional to 5i^q. After the Fourier transformation, D{kyU!lxl'x) = 
e'^'^*' D{kytlxl'x) dt, and some mathematical manipulation we find 



bikytOlX) = D^°HkyLUlX) 

-im{kyu:lx0) 



J 



n -1 



b("\kyLoOQ. (7) 



Here fo is a 2 x 2 identity matrix. It is noteworthy that, except for the matrix form 
and the parametric dependence on ky, Eq. ([7]) is similar in form to the equation 
for Green's function of a crystal with a local defect.^^^ As in this latter equation, 
the first term in the right-hand side of Eq. ([7]) describes bulk excitations, while the 
poles of the second term correspond to excitations localized in the defect region. 

To calculate Green's function _D'^°) {kyUjlxl'^) it is necessary to diagonalize Hamil- 
tonian Q without the term proportional to (5;^o- This can be fulfilled using the 
Bogoliubov-Tyablikov transformation. 



(8) 



where the coefficients Uk^i^k^ and Vkyi^k^ satisfy the usual orthonormality condi- 
tions which follows from the Boson commutation relations of the operators h^yi^ 
and Pkyk^- From these conditions and from the requirement that the Hamiltonian 
be diagonal in the representation of operators (3kyk^ we obtain the following sys- 
tem of equations for the coefficients Ukyi^k^, Vkyi^k^ and the energy of elementary 
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excitations Eh : 



+ -;^[y^ky,l^+l,k^ + V^ky,l^-l,k^ 



(9) 



-Ekyk^v-ky,i,k^ = 2Jv^ky,i^k^ + J cos{ky)ul^i^j. 
J 
2 



with the boundary conditions 

UkyJ^=-l^k^ = 0, V-ky,l^=-l,k^ = 0. (10) 

Sohitions of Eqs. ^ and ^ are standine; waves, 

Kyl^k^ = ^kyk^ S\rL[k^{lx + 1)], 

v-ky,i^k^ = Bkyk^ s\-a{kx{lx + 1)], 

, [2 2J + Ek^k^ 

Akyk^ = \ - " =, (11) 

^ ^ {2J + Ekyk.)^ {2J^kykS 



Bky 



2 2j-fkyk^ 



^ ''2j+EkykS -{^JlkykS 



Ekyk^ — 2 J 

where ^kyk-, ~ \[cos{kx) +cos(fcj,)] and k^ varies continuously in the range (0,7r). 
Using solutions (jlip we find after the Fourier transformation 

D<^'''>{kyU;lX) = f sin[fc,(Z, + 1) sin[fc,(Z; + 1)] 
Jo 

^ ( Pkyk^ Qkyk^ \ ^^^^ 

yuj - Ekyk^ + IT] uj + Ekyk^ + ir] J 

(12) 

p ^ ( ^lyk^ Akyk^Bkyk, \ 

\Akyk^Bkyk^ Bl^^^ J ' 
Qk k — i ^'^y'^' Akyk^Bkyk^ j 

X^Akyk^Bkyk^ A-k^k^ J ' 

where rj = +0. 

The poles of Green's function p2p correspond to bulk excitations - standing spin 
waves (|lip. Apart from them Green's function D[kyUjlxl'^) may have poles connected 
with the second term in the right-hand side of Eq. ([7])- The imaginary parts of 
Green's functions D^^\kyUjlxl'x) and D{kyUjlxl'^) are shown in Figs. [T] and O for 
different distances from the edge. On the edge, Ix — 0, the spectrum lmD{kyUjlxl'x) 



I 



I 
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Fig. 1. The imaginary parts of Green's functions Dn (kyuilxlx) (the solid lines) and D\-^ (kyojlxlx) 
(the dashed lines) for Ix = (a) and Ix = 2 (b). ky = O.Gir. In part (a), the dash-dotted line 
demonstrates the real part of the denominator in the second term in the right-hand side of Eq. JT} • 




Fig. 2. The same as in Fig.[T]for Green's functions Di2{kyLulxlx) and (kyi^lxlx)- 



is dominated by the peak arising from this second term. Indeed, as seen from 
Figs.[lja) and [21^ a), the peak frequency coincides with a zero of the denominator in 



I 



I 
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the term. The standing spin waves manifest themselves as a weak shoulder on the 
high-frequency side of the peak in Fig. [Ija) - the excitations connected with the 
peak eject bulk modes from the edge row of sites. The similar situation is observed in 
the second row. However, as seen in Figs.[ljb) and[2l^b), already in the third row the 
peak is weak and the spectrum is dominated by the continuum of standing waves. 
Thus, the antiferroniagnet is divided into two regions with different spin excitations. 
The two rows near the edge are the location of the mode connected with the pole 
of the second term in Eq. ([7]). This mode has the dispersion uj{ky) = \/2J |sin {ky)\ 
and is termed the boundary spin-wave mode.l^ Near ky = and ky = n it is ill- 
defined - the real part of the denominator of the second term is small but nonzero. 
Excitations of the rest of the crystal are the standing spin waves with dispersion 
(fTT|) . This pictu re is in many respects similar to the situation in the problem of a 
local defect!!^ if local states arise near the defect, they eject bulk states from the 
defect region. 

3. Spin Correlations 

Nearest-neighbor spin correlations can be expressed in terms of correlations of spin- 
wave operators using Eq. Q and the translation invariance of Hamiltonian ([1]) along 
the Y axis, 

(SlSlO = ±^J2{2COS[ky{ly^l'y)]{bk,lJ-k,Mj 

ky 

Mbl,h,Aj + {blrh,,i'j}~\- (13) 

Bearing in mind the property of Green's function ^ Dij{kyUjlxl'x) = Dji{kyUjl'Jx), 
the spin-wave correlations in Eq. (jlSp can be expressed as 

The calculated nearest-neighbor spin correlations parallel and perpendicular to 
the edge, Cl = (Si^+i^j^S;^;^) and Ct = (S/^ j^+iS;^; J, are shown in Fig. g] As 
seen from the figure, main deviations from the bulk value of the correlations fall on 
the edge and the second to the edge row of spins, i.e. on the existence domain of the 
boundary mode. The largest in absolute value spin correlations are observed on the 
edge and between the edge and the second to the edge row. The main contribution 
to these large correlations is made by the boundary mode which has quasi-lD 
character. This result conforms with the fact that the modulus of nearest-neighbor 
spin correlations in the ID antiferromagneJI^ (|(SiSi+i)| = 0.4432) exceeds its 
value in the 2D case (|(SiS]+.a)| = 0.3346 in our calculations, which is close to the 
values obtained earlier P^i^ a is the 2D vector connecting neighbor sites). Thus, 
we relate the enhanced spin correlations on the edge to the separation of the crystal 
into two regions with essentially different excitations - the boundary region with 
the quasi- ID mode and the bulk region with the 2D spin waves. The correlations in 
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Fig. 3. The nearest-neighbor spin correlations parallel (filled circles) and perpendicular (open 
circles) to the edge as functions of the distance Ix to it for T = 0. Bonds perpendicular to the edge 
have their centers at half- integer values of l^- Solid lines are the guide to the eye. The dashed line 
indicates the bulk value. 



the second and between the second and third rows are smaller than in the bulk due 
to the destructive summation of the contributions of these two types of excitations 
[see Eq. d?])]. Qualitatively the obtained picture of spin correlations is similar to that 
observed in the 3D case.L^ However, in the 2D case the deviation of the boundary 
correlations from the bulk value is larger than for 3D - 12% in comparison with 5% 
in the latter case. 

If results in Fig. 3 are compared with the data of Monte Carlo simulations'^ two 
differences stand out: i) in our calculations, the absolute value of the edge correlation 
on the parallel bond is larger than on the bond perpendicular to the edge, while 
in the Monte Carlo data the relation is opposite; ii) there are weak oscillations of 
spin correlations around the bulk value which are perceptible over a few lattice 
periods from the edge in the Monte Carlo data, while such oscillations are missing 
from our results. Partly, these differences can be ascribed to the dissimilarity of the 
used samples - a square-shaped finite crystal in the Monte Carlo simulations and 
a semi-infinite crystal in our case. However, we suppose that the main reason for 
these differences is the quasi-lD character of the boundary mode. This mode is the 
origin of the mentioned peculiarities of spin correlations. However, as known, the 
spin-wave approximation is unsuitable for the ID antiferromagnet. Therefore, it is 
believed that the approach used gives only a rough description of the boundary 
excitations. It would be interesting to compare our results-'' for the 3D case, when 
both the boundary and bulk modes are satisfactorily described by the spin-wave 
approximation, with Monte Carlo data. However, to our knowledge such simulations 
are lacking. 
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4. Conclusion 

In this Letter, we relate the increased spin correlations near the edge of the semi- 
infinite spin-i 2D antiferromagnet to its peculiar spectrum of spin excitations. The 
spectrum involves bulk modes - standing spin waves and the quasi-lD boundary 
mode. This mode is observed in the two boundary rows of sites, and it ejects the bulk 
modes from this region. Thereby the antiferromagnet is divided into two regions 
with different dominant excitations and correlations. The mathematical description 
of the boundary mode and its behavior has much in common with the problem of 
a local defect in a crystal. 
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